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ABSTRACT 

The James Webb Space Telescope (JWST) will enable observations of galaxies at redshifts z > 10 
and hence allow to test our current understanding of structure formation at very early times. Previous 
work has shown that the very first galaxies inside halos with virial temperatures T v i r < 10 4 K and 
masses M vir < 10 8 M at z > 10 are probably too faint, by at least one order of magnitude, to 
be detected even in deep exposures with JWST. The light collected with JWST may therefore be 
dominated by radiation from galaxies inside ten times more massive halos. We use cosmological 
zoomed smoothed particle hydrodynamics simulations to investigate the assembly of such galaxies 
and assess their observability with JWST. We compare two simulations that are identical except for 
the inclusion of non-equilibrium H/D chemistry and radiative cooling by molecular hydrogen. In 
both simulations a large fraction of the halo gas settles in two nested, extended gas disks which 
surround a compact massive gas core. The presence of molecular hydrogen allows the disk gas to 
reach low temperatures and to develop marked spiral structure but does not qualitatively change 
its stability against fragmentation. We post-process the simulated galaxies by combining idealized 
models for star formation with stellar population synthesis models to estimate the luminosities in 
nebular recombination lines as well as in the ultraviolet continuum. We demonstrate that JWST wi\\ 
be able to constrain the nature of the stellar populations in galaxies such as simulated here based on 
the detection of the Hel640 recombination line. Extrapolation of our results to halos with masses both 
lower and higher than those simulated shows that JWST may find up to a thousand star-bursting 
galaxies in future deep exposures of the z > 10 universe. 

Subject headings: cosmology: observations - galaxies: formation - galaxies: high-redshift - hydrody- 
namics - intergalactic medium - stars: formation 



1. INTRODUCTION 

The hierarchical assembly of dark matter halos and 
the cooling and condensation of the cosmic gas to form 
stars and galaxies inside them dRees fc Ostriker 119771 : 
ISilkl [19771 : IWhite fc Reelfl978l : iBlumenthal et al.l[l98l 
are major pillars of the current cold dark matter (CDM) 
paradigm of structure formation in the universe with 
cosmological consta nt A. Both (s e mi-) analytical ar- 
guments (e.g.. iTegmark et al.l 119971 : iReed et all 120051: 
iNaoz et"al1 120061) and simulations (e-g-, IBromm et all 
120021 : lAbel et al.ir2QQ2t lYoshida et aTll2006l ) suggest that 
the first stars have formed at redshifts as high as z ~ 
30, when the universe was just abo ut a percent of 
■ its present age (for reviews see, e.g., Barkan a" fc Loebl 
120011 : IBromm fc Larsonll2004t IBromm et al.l2QQ9h . Their 
light ended the Dark Ages that followed the release 
of the cosmic microwave background (CMB) radia- 
tion at z ~ 1100 and fundamentally transformed the 
universe during a landmark period called the epoch 
of re io nization (for reviews se e , e.g., lLoe b & Ba rkanal 
2001 ICiardi fc Ferraral [2005[ iB arkana fc Loebl 120071 : 
Trac fc Gnedinll2009l : IStiavellill2QQ9l : lLoebll2010[ r^ 

Future observations with telescopes such as, for ex- 
ample, Planck 1 , the Low Frequency Array 2 , the Murchi- 
son Widefield Array 3 , the Atacama Large Millime- 
ter Array 4 , and the James Webb Space Telescope 

1 sci.esa.int/planck/ 

2 http://www.lofar.org 

3 http:/ /www. haystack. mit.edu/ast/arrays/mwa/ 

4 http:/ /www. almaobservatory.org/ 



(JWST) 5 will test our current theoretical understand- 
ing of the formation of stars and galaxies at these 
early times. A fascinating prospect is the detection 
of line and continuum radiation from the first galax- 
ies with JWST. The ratio of the hydrogen and helium 
Balmer line luminosities from recombining gas has been 
proposed as a telltale signature that distinguishes be- 
tween first-generation, metal-free (Population III) and 
subsequent metal-enriched (Population II) star forma- 
tion or between stellar s ources and b l ack holes (e.g., 
Tumlinson fc Shull 2000; Bro mm et all [200 ll lOh et al.l 
2001 1 : iTumlinson et al.l 120011 : ISchaererl [20021 : iSchaererl 
2003; Uohnson et al.ll2009f ). In addition, the detection of 
Lya, molecular, or metal line cooling radiation from high 
redshifts would probe the gravitational assembly of the 
gas in the first halos (e.g., lHaiman et al. 2000; iDijkstral 
2009, Mizusawa et al. 2005; Appleton et al. 2009),' offer- 
ing direct insights in the structure and dynamics of the 
first galaxies and the surrounding intergalactic medium 
(e.g., Santos 2004; Diikstra et al. 2006; Verhamme et al. 
l2006UDiikstra et al.l 120071 : lLaursen et al.ll2010D . 

The combination of upcoming observations with JWST 
and other future telescopes with detailed numerical su- 
percomputer simulations of the first galaxies and reion- 
ization will transform our knowledge of structure for- 
mation in the universe. Most of the numerical effort 
in early galaxy formation has concentrated on investi- 
gating the properties of the very first building blocks 
of galaxy assembly, minihalos and dwarf galaxies with 

5 http://www.jwst.nasa.gov/ 
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virial temperatures < 10 4 K, 
masses < 1Q 8 at z > 10 (e.g 
Bromm et all l2QQ2b iWise et al.l l2Qq 



corre sponding to halo 



Abel et al. 


2002 


Turk et al. 


2009; 



TABLE 1 
Simulation Parameters 



Stacv et al.ll2Q10b iGreif et al.ll2010D .~A key result emerg- 



ing from the existing work is that the luminosities of 
these very low-mass objects are, unless magnified by 
gravitational lensing, too low, by at least one order of 
magnitude, t o be detected in e ven very deep expos ures 
with JW ST (IGreif et al.l [20091: I Johnson et al.l [20091: 



also, e.g. , lLoeb & Haiman 1997; Hai man fc Loebl 119981 : 
|pail99a lOh et al.ll200ll : iTrenti et alJl2QQ~) . The light 
collected in future deep-field observations with the JWST 
may thus be dominated by emission from dwarf galax- 
ies inside halos that are about ten times more massive, 
M vir - 10 9 M . 

Our goal is to extend and complement existing nu- 
merical work on the first galaxies by investigating the 
role played by galaxies inside halos with masses M v i r > 
10 9 Mq, at times before and during the epoch of reion- 
ization, when these galaxies were assembling, possibly 
contr i butin g a significant fract i on (if not most ; e.g. , 
lLoebl 120091: IWise fc CerJ 120091 : iSalvaterra et al.l 120101 : 
iRaicevic et al.ll201Q[ ) of the ionizing emissivity in the uni- 
verse, to the present day, when these galaxies may be 
found in the Local Group as fossil probes of the begin- 
nings of galaxy formation (for reviews see, e.g., iMateol 
[1991 iTolstoy et al.ll200l iRicottil 1201(1 lMaverll20loh . In- 
deed, the new field of 'dwarf archaeology' may hold the 
key to unravel the interplay of star and galaxy formation 
at th e end of the cosmic dark ages (Fre bel fc Bromml 
|201Q[ ). Here, we report our first steps towards achieving 
this goal by studying the assembly of dwarf galaxies in 
halos reaching virial masses M vir - 10 9 M at z = 10 us- 
ing cosmological smoothed particle hydrodynamics sim- 
ulations. 

We utilize a zoomed simulation technique that allows 
us to simulate the gravitational and hydrodynamical pro- 
cesses of dwarf galaxy formation at high resolution while 
keeping information about structure formation at large 
representative scales. Our simulations include radiative 
cooling from atoms and molecules in gas of primordial 
composition but ignore star formation and the associated 
feedback. We post-process our simulations with idealized 
models for star formation and employ population synthe- 
sis models to estimate the prospects for a direct detection 
of the first galaxies with the upcoming JWST. Other as- 
pects of the simulated galaxies, like e.g., their role as 
reionization sources, will be investigated in subsequent 
works, in which we will explicitly account for the effects 
of star formation and associated feedback. 

The structure of this paper is as follows. In Section [2] 
we describe the set-up of our simulations. Then, in Sec- 
tion [3j we present our results, and subsequently discuss 
the assembly of the simulated halo, its structure and dy- 
namics at redshift z = 10 and the properties of the disks 
it hosts. In Section [4] we combine our simulations with 
assumptions about star formation to assess the observ- 
ability of the first galaxies in future observations with 
JWST. In Section [5] and Section [6] we discuss, respec- 
tively, implications and limitations of our work. Finally, 
in Section [71 we summarize our work. 

Throughout this work we assume ACDM cosmological 
parameters ^ m = 0.258, ^ b = 0.0441, ^ A = 0.742, a 8 = 
0.796, n s = 0.963, and h = 0.719, which are con- 
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sisten t with the 5-year (jKomatsu et al.l l2QQ9h and 7- 
year (Komat su et al.ll2010l ) analyses of the observations 
with the Wilkinson Microwave Anisotropy Probe satel- 
lite. Distances are expressed in physical (i.e., not co- 
moving) units, unless noted otherwise. 

2. SIMULATIONS 

We use a modified version of the N-body/TreePM 
Smoothed Part i cle Hydrodynamics ( SPH) code gadget 
([Spri ngel 2005|; iSpringel et al.l l200ll ) to perform a suite 
of high-resolution zoomed cosmological hydrodynamical 
simulations in a box of size L = 3.125 h~ x Mpc comoving. 
The box size was chosen after inspec tion of published 
dark matter halo mass functions (e.g., Reed et al. 2007) 
such that the box contains at least one halo of mass 
- 10 9 M at redshift z = 10. 

We carry out SPH simulations of primordial metal- 
free gas including non-equilibrium radiative cooling from 
both molecular and atomic species and from atomic 
species only. In the following, these two types of simu- 
lations will be distinguished by an additional NOMOL 
at the end of the label of the atomic cooling simula- 
tion. The simulations are performed with the grav- 
itational forces softened over a sphere of Plummer- 
equivalent radius e. Our simulations Z4 and Z^NOMOL, 
which are obtained by zooming into a parent cosmo- 
logical simulation 4 times, use a force softening radius 
e = 0.1 h~ x kpc comoving applied to all particles. 

We employ the entropy- co nserving formulation of SPH 
(jSpringel fc Hernquistll2002[ ) with TVneigh = 48 neighbor 
particles per SPH kernel. We limit the radius h of the 
SPH kernel to above a fraction fh of the softening length, 
h > fht, where fh = 0.01. The simulations are summa- 
rized in Table [TJ 

2.1. Initial Conditions 

All simulations start at redshift z = 127. Initial parti- 
cle positions and velocities are obtained by applying the 
Zel'dovich approximation ([Zerdovichl [1970) to particles 
arranged on a Cartesian grid. We adopt a transfer func- 
tion for matt er perturbations generate d with cmbfast 
(version 4.1; iSeljak fc Zaldarriagalll996l ). 

To achieve high resol ution we make use of the zoomed 
simulation technique (Navarro fc White! 11995 we use 
the same code as in Grei f et al.l 12008). We first per- 
form a simulation in which the initial conditions are set 
up using 2 x 64 3 (dark matter and gas) particles ar- 
ranged on a uniform Cartesian grid. We then use the 
friends-of-friends (FOF) halo finder, with linking param- 
eter b = 0.2, built int o the substructure finder subfind 
([Springel et al.ll2001bl ) to locate an FOF halo with mass 
> 10 9 M at z = 10 in this simulation. Using subfind, 
we determine the most bound particle of this FOF halo 
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Fig. 1. — Halo assembly history in the ZJ± simulation. Top panel: 
Growth of the virial mass M v j r associated with the most massive 
progenitor FOF halo (solid curve) is close to exponential (dotted 
line). Bottom panel: Accretion rate corresponding to the halo 
growth shown in the top panel. The accretion rate is consistent 
with an alytical estimates of the mass growth rate of the main pro- 
genitor (Neistein et al. 2006; dash-dotted curve). 

and let it mark the center of the region that we wish to 
resimulate. We compute the associated virial radius r v i r 
by determining the radius of the sphere around the most 
bound particle within which the average matter density 
equals 200 times the critical density at z = 10. The par- 
ticles within a region of radius r < 3r v i r around the most 
bound particle are traced back to their locations in the 
initial grid where they mark the region of refinement. 

Subsequently, all parent particles in a cube enclosing 
the refinement region are replaced by 8^ daughter par- 
ticles, where N\ is the zoom level. Our simulations Z4 
and Z4NOMOL employ N\ = 4 and hence in these sim- 
ulations the daughter gas (dark matter) particles have 
masses m g = 484 M (raDM = 2350 M©). To reduce 
numerical artifacts due to the large difference in the par- 
ticle masses for particles inside and outside the refine- 
ment region (mass ratios 8^), the refinement region is 
surrounded by N\ — 1 concentric, nested layers in which 
the parent particles are successively replaced by S^ 1-1 , 
S^ 1-2 , 8 daughter particles and, hence, the particle 
masses vary gradually, by discrete factors of 8, with in- 
creasing distance to the refinement region. The simula- 
tion is then re-run after applying the Zel'dovich approx- 
imation to evolve all particles to the starting redshift 
z = 127. 

2.2. Chemistry and Cooling 

All our simulations include radiative cooling in the op- 
tically thin limit. We assume that the gas is of primordial 
composition using a hydrogen mass fraction X = 0.752 
and a helium mass fraction Y = 1 — X. We follow the 



non-equilibrium chemistry and cooling of H2, D, HD, D + , 
H + , H, D, and He, and we incl ude H and Hj~ assuming 
their equilibrium abundances (j Johnson fc Bromml l2006: 
IGreif et al.ll20loh . 

In simulation Z4, gas cools by collisional ionization and 
excitation, the emission of free-free and recombination 
radiation, Compton cooling off the CMB, and emission of 
radiation by molecular hydrogen and hydrogen deuteride 
(HD). If initial abundances are expressed as number den- 
sity with respect to hydrogen, where nn = Xp g /mn 
with p g being the gas density at z = 127, and ran is 
the mass of the proton, we choose: H2 = 1.1 x 10 -6 , 
D = 2.6 x 10- 5 , HD = 10- 9 , D+ = 1.2 x 10~ 8 , 
H+ = 3 x 10- 4 , He+ = and He++ = 0, from which 
the initial abundances of the remaining species (H, D, 
He) as well as the abundance of electrons are obtained 
through application of conservation laws. Our choices for 
the initial abundances are consistent with computations 
of cosmological abundances in the ear ly universe (e.g., 
iLepp fc Shuflir^lGamfc Pallalll998h . 

Thanks to molecular cooling, gas in simulation Z4 
may reach temperatures as low as ~ 10 2 K. Sim- 
ulation Z4NOMOL is identical to simulation Z4 ex- 
cept that the formation of molecular hydrogen is sup- 
pressed, as would be the case in the presence of a strong 
photo-dissociating Lyman- Werner radiation ba ckground 
(jStecher fc Williams! [19671 : lHaiman et all 119971 ). Gas in 
simulation Z4NOMOL therefore cools only via atomic 
processes, which are inefficient in reducing the ther- 
mal energy of primordial gas with temperatures below 
- 10 4 K. 

2.3. Jeans Floor 

Simulations with mass resolutions insufficient to re- 
solve the Jeans mass Mj = (47r/3)p m (Aj/2) 3 by at 
least N res = 2 SPH kernel masses Mk = ^neigh^e- 
may suffer from artificial fragmentation (jBate fc Burkertl 
1 19971 ). Here, p m is the total (dark matter and gas) 
mass density, Aj = c s 7r 1 / 2 (Gp m ) -1 / 2 the Jeans length, 
c s = [^yk^T / {firriYi)} 1 ^ 2 the adiabatic speed of sound, 7 
the ratio of specific heats, and \i the mean gas particle 
mass in units of the proton mass. Our finite mass reso- 
lution implies a maximum density 

n H , ma x = 8 x 10 s cm- 3 / g (£) (^) (1) 

up to which we satisfy the IB ate fc Burkertl ([19971 ) crite- 
rion, where f g = p s /p m - 

To satisfy the (Bate fc Burkertl (|1997l ) criterion inde- 
pendent of densit y we make use of a density-d ependent 
temp erature floor (Robertson & Kravtsov 2008; see also, 
e.g., ISchave fc Da lla Vecchia 20081 for a related ap- 
proach). At each time step and for all gas particles we 
compute the Jeans mass Mj and compare it to the reso- 
lution mass N tqs Mk. If the Jeans mass becomes smaller 
than the resolution mass, then we increase the parti- 
cle internal energy and, hence, the particle temperature 
such that the Jeans mass becomes equal to the resolution 
mass. 

3. RESULTS 

In this section we describe the outcome of our simula- 
tions. We start in Section 13.11 by briefly presenting the 
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Fig. 2. — Mass and density profiles centered on the most bound halo particle in simulation ZJ^. In all panels, the vertical lines on the left 
mark the gravitational softening radius e and the vertical lines on the right mark the virial radius at z = 10, the final simulation redshift. 
Left panel: Dark matter (red dash-dotted curve) and gas density profiles at z = 10. The gas densities were computed both by dividing 
the gas mass summed in shells by the shell volumes (blue dashed curve; this was also done for the dark matter density profiles) and by 
averaging the SPH particle densities inside shells (green dotted curve). For comparison, the corresponding SPH particle density profile from 
simulation Z4NOMOL is also shown (black dotted curve). The dark matter profile is approximately singular-isothermal for all distances 
larger than the softening radius. The two steepenings in the gas density profiles at r « 0.07 kpc and r « 0.3 kpc are due to the presence of 
two nested disks. Middle panel: Enclosed gas, dark matter and total masses at z = 10. Right panel: Evolution of the enclosed gas mass. 
There is a rapid inflow of mass into the central unresolved region r<eatll.5<z< 12.5. 



growth histories of the simulated halos. We focus our 
subsequent dis cussion on the halo properties at z = 10. 
In Section 13.21 we discuss the structure of the halos, in 
Section 13.31 we describe the dynamics of the gas inside 
them, and in Section 13.41 we investigate the emergence 
of nested gas disks at the halo centers. Throughout we 
will discuss differences and similarities between simula- 
tion Z4 and simulation Z4NOMOL in which molecular 
hydrogen formation is suppressed. 

3.1. Growth 

We use FOF to locate the simulated halos at the final 
simulation redshift z = 10, and then compute the halo 
properties as follows. Given a FOF halo, we use subfind 
to identify its most bound particle and let it mark the 
halo center. We then compute the virial radius of the 
halo by determining the radius of the spherical volume 
centered on the most bound particle within which the 
average matter density is equal to 200 times the critical 
density at z = 10. The total mass inside this volume 
defines the halo virial mass. 

We find that r vir « 3.1 kpc and M vir w 1.3 x 10 9 M , 
independent of the inclusion of molecular cooling. For 
the adopted cosmological parameters, this virial mass 
corresponds to « 3a fluctuations in the linear theory 
density field (e.g., Barka na~fc Loebl [2001). The circu- 
lar velocity v c = (GM^/r^) 1 / 2 and virial temperature 
T v i r = /imn^vir/ (3&b) implied by the virial mass and the 
virial radius are v c ~ 40 km s _1 and T v i r « 42000 K, 
where we have assumed /i = 0.6 appropriate for ionized 
gas with primordial composition. 

After having located the halo at z = 10, we trace its 
history to higher redshifts using the FOF halo finder to- 
gether with subfind, as follows. Knowing the FOF halo, 
the descendant, at redshift Z{ corresponding to simu- 
lation snapshot z, we locate the FOF halo at redshift 
Zi-i > Z{ corresponding to snapshot i — 1 that shares, 
among all FOF halos present at 2^_i, the most mass 
with the descendant. We then use subfind to identify 
the most bound particle within this FOF halo and ob- 



tain the properties of the halo at Zi-\ by computing its 
virial radius and mass in a sphere of average matter den- 
sity 200 times the critical density at 2^_i centered on the 
most bound particle. 

Figure [T] shows the evolution of the halo mass M v i r (top 
panel) and its corresponding rate of growth dM w \ v /dt 
(bottom panel) in simulation Z4. The mass growth 
is described well by an exp onential fit M v { r (z) = 
M exp(-az) (jWechsler et all l2002f ) with M = 2 x 
10 11 M and a = 0.5. The derived growth rates are 
consistent with analytic al estimates of the ra te of growth 
of the main progenitor (jLacev fc C ole 1993). The dot- 
dashed c urve shows the growt h rate as given in equation 
(A15) of iNeistein et al.l (|2006f ) with q = 2.3. 

3.2. Structure 

The left panel of Figure [2] shows density profiles spher- 
ically averaged around the most bound halo particle ob- 
tained from simulation Z4 at z = 10. The dark matter 
density profile (red dash-dotted curve; obtained by sum- 
ming particle masses inside spherical shells and dividing 
by the shell volume) follows an isothermal shape p oc r~ 2 
for radii e < r < r v i r . The gas density profile does not 
follow the shape of the dark matter profile but shows sig- 
nificant small-scale structure. We have computed the gas 
density profile both by averaging the SPH particle den- 
sities inside spherical shells (green dotted curve) and by 
summing particle masses inside spherical shells and di- 
viding by the shell volume (blue dashed curve). The two 
methods for computing the density profile yield different 
results because the spatial distribution of the gas mass is 
highly non-uniform, as will be discussed below. We also 
show the gas density profile in simulation Z4NOMOL 
(black dotted curve). The dark matter density profile 
obtained in simulation Z4NOMOL is nearly identical to 
that from simulation Z4 and hence is not shown. 

The middle panel of Figure [2] shows the enclosed mass 
as a function of distance r from the halo center for sim- 
ulation Z4. The cumulative mass of dark matter dom- 
inates the significantly more centrally concentrated cu- 
mulative mass of gas for radii r > 0.4 kpc. The build-up 



THE FIRST DISK GALAXIES 



5 




Fig. 3. — Densities (top) and temperatures (bottom) at z = 10 in simulation Z4 in cubical slices centered on the most bound halo particle. 
The left panels present edge-on views of the disks and encompass a volume slightly larger than the virial region with radius r v i r ~ 3.1 kpc. 
The middle panels are zooms into the cubical regions marked in the left panels. The right panels are zooms into the cubical region marked 
in the middle panels and are reoriented such that the outer disk is seen face on. The temperature of the underresolved gas is artificially 
elevated because of the use of a density-dependent temperature floor to prevent artificial fragmentation. The spirals do not show signs of 
fragmentation; the gas clump seen in the bottom of the face-on view of the disk (right panels) is a gas-rich subhalo in projection. 
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Fig. 4. — Same as Figure [3] but for simulation Z4NOMOL. As in Figure [3] the left-hand and middle panels present edge-on views of the 
galaxy, and the right-hand panels show the galaxy face-on. Because of the lack of efficient low-temperature coolants and in contrast to 
simulation Z4 (Figure [3}, the gas inside the filaments is at about the same temperatures as the diffuse gas in between the filaments. The 
temperature of the underresolved gas core is artificially elevated because of the use of a density-dependent temperature floor to prevent 
artificial fragmentation. 
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Fig. 5. — Molecular hydrogen fraction ?7h 2 = nn 2 /nn at 2 = 10 
in simulation ZJ^.. The image shows the same cubical region as the 
left panels in Figure [3] The dashed circle marks the virial radius 
r v i r ^3.1 kpc. The gas inside filaments and subhalos has large 
molecular fractions ?7h 2 ~ 10 -3-5 — 10 -3 which enables efficient 
low-temperature radiative cooling. 



of the final gas mass profile shown in the middle panel 
is illustrated in the right panel of Figure [2] There is 
a large rapid increase in the mass of the central region 
(r < e) around redshift z & 12. In little more than 
40 Myr (11.5 < z < 12.5), the central gas mass grows 
from - 2 x 10 7 M to - 5 x 10 7 M . In Section g] we 
will assume that this rapid collapse of large gas masses 
triggers a massive burst of star formation in the cen- 
tral core. The cumulative mass profiles from simulation 
Z4NOMOL exhibit a nearly identical behavior and again 
are not shown. 

Figures [3] and |4] give further impressions of the bary- 
onic structure of the simulated halos at the final simu- 
lation time, i.e., at redshift z = 10. All quantities are 
shown for both simulation Z4 (Figure [3]) and simulation 
Z4NOMOL in which molecular hydrogen formation was 
suppressed (Figure Hj). The panels in the left columns 
of Figures [3] and |4] show the hydrogen number densities 
and temperatures, mapped to a three-dimensional grid 
using standard mass-conserving SPH interpolation and 
mass- weighted averaging along the line of sight, within a 
cubical volume encompassing the virial region. 

The panels show that independent of the inclusion 
of molecular cooling, gas is organized in four geomet- 
rically distinct components: diffuse low-density (nn < 
10 -2 cm -3 ) gas, collimated streams, or filaments, of 
smooth dense (rtn < 10 _1 cm -3 ) gas that penetrate deep 
into the virialized region, dense (nn > 10 _1 cm -3 ) gas- 
rich clumps with mostly-spherical appearance and a cen- 
tral dense (nn > 10 1 cm -3 ) gaseous flattened object seen 
edge-on. The dense clumps inside the virial radius are 
associated with low-mass (< 10 7 M ) halos that entered 
the virial region before z = 10. In the following we refer 
to these halos as subhalos. 

The panels in the middle columns of Figures [3] and |4] 
are zooms into the cubical regions marked by the white 
solid rectangles in the panels of the left columns. The 
panels in the right columns are zooms into the cubical 
region marked in the panels of the middle columns, but 
with their coordinate axes rotated. The zooms resolve 



the central flattened object into two nested disks whose 
orientations are tilted with respect to each other. The 
disks cause the steepenings of the gas density profile at 
r « 0.07 kpc and r « 0.3 kpc seen in the left panel of 
Figure [21 We let these radii define the sizes of the disks. 
Note that the disks surround a central unresolved core of 
radius r < e. The disks will be discussed in more detail 
in Section [3T4l 

The images of the gas temperatures reveal a qual- 
itative difference between simulation Z4 and simula- 
tion Z4NOMOL in which molecular hydrogen formation 
was suppressed. While in both simulations the ten- 
uous gas in between the filaments is at temperatures 
close to the virial temperature of the halo, the tem- 
peratures of the gas inside filaments, subhalos and the 
disks are up to an order of magnitude lower in Z4 than 
in Z4NOMOL. Figure [5] shows the molecular hydrogen 
fraction t]yl 2 = n H 2 / n H inside the same cubical region as 
shown in the left panels of Figure 03 The molecular frac- 
tion is greatly increased up to T]u 2 ~ 10 -3 in gas with 
densities > 1 cm -3 and temperatures T < 1000 K, 
which resides mostly in the filaments and subhalos and 
in the central halo region. 

Without sufficient molecular hydrogen, the intra-halo 
gas cannot cool efficiently below temperatures T ~ 10 4 K 
as atomic cooling is exponentially suppressed because of 
the lack of thermal excitation of bound electrons. Gas 
that is shock- heated to T > 10 4 K upon entry in the 
halo then remains hot at T w 10 4 K until it is incorpo- 
rated into the dense disks where it cools to slightly lower 
temperatures. In the presence of a sufficient amount of 
molecular hydrogen, on the other hand, radiative cooling 
counters virial heating in the dense gas inside filaments 
down to temperatures T < 10 3 K. The ac cretion along 
the filaments then occurs in a cold mode (Wise & Abell 
20071: iGreifet al]l2QQ8t see also iBirnboim feDekell 120031: 
Keres et al.ll2005l : iBrooks et al.ll2009l : Ivan de Voort et al.l 
20101 for cold accretion of gas inside more massive halos). 



3.3. Dynamics 

Figure [6] shows spherically averaged (mass- weighted) 
profiles of the gas radial velocities v gjY (left panel) and 
fractional radial velocities v gjY /v g (right panel) at z = 10 
in simulation Z4, where v g = |v g | and v g is the gas ve- 
locity. For comparison, the corresponding velocities for 
the dark matter are also shown. The velocities were cor- 
rected for the bulk halo motion by subtracting the veloc- 
ity of the center of mass of all gas particles within the 
virial region and were calculated relative to the location 
of the most bound particle. Figure [6] shows that both 
the dark matter and the gas approach the virial region 
along mostly radial orbits with similar infall velocities 
consistent with the halo circular velocity (Section 13. 2p . 
Inside it, their velocity distributions, however, differ sig- 
nificantly. 

The dark matter isotropizes just upon entry in the 
virial region, i.e. at r « r v i r (vertical line on the right), 
which is reflected in a sharp drop of the ratio of radial to 
total velocity to 3 -1 / 2 expected for an isotropic velocity 
distribution (horizontal dotted line in the middle panel 
of Figure [6|). The gas, being able to radiatively cool and 
lose gravitational energy, keeps streaming with radial ve- 
locities that increase towards the halo center and reach 
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Fig. 6. — Gas flow at z = 10 around the most bound halo particle in simulation Z4- In all panels, the vertical line on the left shows 
the gravitational softening radius e and the vertical line on the right marks the virial radius r v [ Y . Left panel: Spherically averaged radial 
particle velocities for dark matter (blue dashed curve) and gas (black solid curve). The horizontal line marks zero radial velocity. Middle 
panel: Same as left panel, but with particle radial velocities divided by the particle total velocities. The horizontal line marks the ratio of 
radial and total velocity expected for an isotropic velocity distribution. While the dark matter isotropizes after entering the virial region, 
the gas keeps falling in along mostly radial orbits until it reaches the outer disk at r < 0.3 kpc. Right panel: Gas accretion rates. For 
comparison, we also show the gas accretion rates in simulation ZJ^MOMOL. The spikes are associated with gas-rich subhalos, a majority of 
which accretes along filaments. Note that the scale of the vertical axis is linear. 
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Fig. 7. — Radial (top row) and line-of-sight (i.e., along the axis of projection; bottom row) gas velocities at z = 10 in simulation Z4- 
The images correspond in size and orientation to the images in Figure [3] which allows the matching of features in velocity and real space. 
As in Figure [3] the left-hand and middle panels present edge-on views of the simulated galaxy, and the right-hand panels show the galaxy 
face-on. The top row panels show that the central disks grow through infall of dilute gas, channeling of both dense smooth gas and gas-rich 
subhalos along filaments and merging with gas-rich subhalos from outside filaments. The bottom row panels show the kinematic signature 
of the disks. 
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maximum values — v gyr < 60 km s _ . The gas eventually 
hits the central disk at r « 0.3 kpc where it circularizes. 
The right panel in Figure [6] shows that at redshift z = 10, 
in both simulation Z4 and simulation Z4NOMOL, the 
spherically integrated gas accretion rates are, on aver- 
age, — r 2 f p g v g ^ r dft « 1 M yr _1 , independent of radius 
r > 0.07 kpc. The prominent spikes exhibited by the gas 
accretion rates are due to the infall of gas-rich subhalos 
along the radial filaments. 

The complexity of the gas dynamics within the virial 
region in simulation Z4 is revealed by the radial and line- 
of-sight (i.e., along the axis of projection) velocities that 
are shown, respectively, in the top and bottom panels 
of Figure [71 The velocity structure of the gas in sim- 
ulation Z4NOMOL is very similar. The images show 
views that correspond in size and orientation to the views 
shown in Figure [3l This allows to identify structures in 
the velocity images with those in the images of density 
and temperature. The velocity images were obtained by 
mapping particle radial and line-of-sight velocities to a 
three-dimensional grid using standard SPH interpolation 
and performing a mass-weighted average along the line 
of sight to project the grid into a two-dimensional plane. 

The line-of-sight velocities (bottom panels in Figure [7]) 
show the kinematic signatures of two rotating disks that 
are centered on a spatially unresolved core. The inner 
disk continues to grow in mass not only through accre- 
tion of gas from outside the disk plane but also through 
accretion of gas from within the outer disk to which it 
is kinematically connected (bottom middle panel in Fig- 
ure [7j). The radial velocity (top panels in Figure [7]) shows 
a complex inflow pattern. The gas stream that points to 
the top right corner in the top middle panel of Figure [71 
is gas ejected by the passage of a gas-rich subhalo from 
beneath the disk immediately before z = 10. A major- 
ity of subhalos appears to be associated with filaments, 
which resembles the picture of anisotropic accretion of 
satellites in simulatio ns of massive galaxies and clusters 
at low redshift (e.g., iLibeskind et al.l [2010b iKnebe et al.l 
12001 . 

3.4. The Disks 

Figure [8] shows the assembly of the inner and the outer 
disk in simulation Z4- The disk assembly times and histo- 
ries in simulation Z4NOMOL are very similar. The inner 
disk forms at z ~ 13.5 as a result of a major merger. Ini- 
tially, it is relatively thin and extended and also shows 
marked spiral structure. Its gas, however, collapses into 
the highly concentrated disk that is seen at z = 10 in Fig- 
ure [3] rather quickly, within a redshift interval Az < 1, 
corresponding to < 40 Myr. The outer disk is assembled 
from the halo gas at z ~ 11.5, i.e., < 100 Myr after the 
assembly of the inner disk. It grows in size and develops 
increasingly pronounced spiral structure. In none of our 
simulations do the disks show signs of fragmentation. 

The right panels in Figures [3] and [4] present face-on 
views of the disks at z = 10. In simulation Z4, the outer 
disk shows much more developed spiral structure. When 
seen edge-on (middle panels of Figures [3] and [4]), the 
outer disk in Z4 looks somewhat thinner and more per- 
turbed than the outer disk in Z4NOMOL. The latter may 
be explained in part by the fact that in simulation Z4, 
thanks to the efficient low-temperature molecular cool- 
ing, subhalos have larger gas fractions than in simulation 



Z4NOMOL, which increases the chance for possibly vio- 
lent subhalo-disk interactions. 

In simulation Z4NOMOL, the disk gas reaches mini- 
mum temperatures T > 8000 K slightly below the tem- 
peratures in the diffuse gas and filaments because the in- 
creased densities in the disks imply shorter cooling times. 
In simulation Z4, on the other hand, the disk gas can 
cool to temperatures T < 1000 K thanks to the pres- 
ence of molecular hydrogen. Note that the disk tem- 
peratures in simulation Z4 are also determined by the 
temperature floor enforced to prevent artificial fragmen- 
tation (Section |2.3|) . The temperature floor affects the 
evolution of the gas for densities above nn > 10 cm -3 
in simulation Z4 and densities above nn > 10 6 cm -3 in 
simulation Z4NOMOL (see equation pQ). These densi- 
ties correspond, respectively, to radii r < 0.3 kpc and 
r < 0.03 kpc (see Figure [2]). 

The final mass distributions in the disk region in the 
simulations Z4 and Z4NOMOL are very similar. In both 
simulations, the volume inside r < 0.3 kpc contains gas 
and total masses of « 1.4 x 10 8 M and « 2.5 x 10 8 M , 
respectively (see Figure [2]). These masses amount to ~ 
65% of the gas mass and to ~ 25% of the total mass 
inside the virial region. In simulation Z^, roughly 20% 
(35%) of the total (gas) mass within r < 0.3 kpc is in 
the central unresolved core, about 30% (40%) is in the 
inner disks and about 50% (25%) is in the outer disks. In 
simulation Z4NOMOL, roughly < 17% (< 28%) of the 
total (gas) mass out to radii r < 0.3 kpc is in the central 
unresolved core, < 32% (< 40%) is in the inner disk and 
about 51% (< 32%) is in the outer disk. 

Figure [9] shows several azimuthally averaged properties 
of the disks at z = 10 in simulation Z4 and Z4NOMOL. 
The left panel of Figure [9] shows rotational gas veloc- 
ities Vg 5r ot (black curves), adiabatic sound speeds c s = 
[jk B T/(fim u )} 1/2 (blue curves) and radial velocities v v r 
(red curves), all with respect to the Keplerian veloc- 
ities vk = [GM(< r)/r} 1 / 2 and for both simulation 
Z4 (solid curves) and simulation Z4NOMOL (dashed 
curves). The rotational velocities were computed using 
^g,rot = ( v g — Vg,!-) 1 ^ 2 - For computing the sound speed 
we have assumed a ratio of specific heats of 7 = 5/3 
appropriate for a mostly atomic gas and assumed that 
the cold disk gas is mostly neutral, i.e., \i — 1.2. Both 
the inner and the outer disks exhibit a high degree of 
rotational support with the outer disk showing nearly 
Keplerian motion. The rotational velocities are larger 
than the sound velocities by factors of < 5 — 10. Radial 
velocities are small compared to rotational velocities at 
all disk radii and drop sharply to zero once the gas hits 
the central core at r < e. 

The middle panel of Figure [9] shows the gas surface den- 
sity profiles. In both Z4 and Z4NOMOL, the inner disk, 
which is resolved with > 5 gravitational softening radii, 
is characterized by an exponential surface density pro- 
file with scale length 0.015 kpc that extends over several 
scale lengths. In Z4NOMOL, the outer disk follows an 
exponential profile over several scale lengths of 0.05 kpc. 
In Z4, the surface density profile shows significant devia- 
tions from an exponential profile. The nearly exponential 
density profiles of the outer disks in our simulations are 
only gradually built up and at higher redshifts the den- 
sity profile in simulation Z4 shows pronounced ripples 
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Fig. 8. — Assembly of the inner and the outer gas disks in simulation Z4- The panels present face-on views of the simulated galaxy and 
show the evolution of gas densities in the redshift range 11.0 < z < 14.5 in cubical slices of linear size 0.5 kpc. They can be compared to 
the top right panel of Figure [3] which shows a face-on view of the gas density at the final simulation redshift in the same cubical slice. The 
color coding is identical to that in Figure [3] The outer disk forms at z ~ 11.5, roughly 100 Myr after the assembly of the inner disk at 
z % 13.5. 




Fig. 9. — Properties of the disks at z = 10. Left panel: Rotational gas velocities v g , r ot (black curves), adiabatic sound speeds c s (blue 
curves) and radial gas velocities v g , r (red curves) with respect to Keplerian velocities in simulations Z4 (solid curves) and Z4NOMOL 
(dashed curves). Middle panel: Gas surface density profiles in simulations Z4 (black solid curve) and Z4NOMOL (blue dashed curve). 
Right panel: Toomre Q parameter in simulations Z4 (black solid curve) and Z4NOMOL (blue dashed curve). In each panel, the vertical 
line marks the gravitational softening radius e. The gas in the outer disk (0.1 kpc < r < 0.3 kpc) is on nearly Keplerian orbits. Both the 
inner (e < r < 0.07 kpc) and the outer disk show surface density profiles that are roughly exponential with scale lengths as indicated in 
the middle panel (dotted lines). Toomre parameters Q > 1 indicate stable disks. 



due to the presence of thick spiral structure. 

The right panel of Figure [9] shows the Toomre pa- 
rameter Q = c s ft/(7rGT), where c s is the velocity 
dispersion, n = (4^ 2 + rdVt 2 / dr) 1 ^ 2 is the epicyclic 
frequency and Vt = v g:rot /r is the angular velocity 
([Toomrd 19641). A st a ndard linear theory instability 
analysis (IToomre 1964; Goldreich Lynden-Bell 119651 : 
iBinney fc Tre maine 2008) shows that a gas disk becomes 
unstable to axisymmetric perturbations if Q < 1; the pre- 
cise threshold for instability depends on the disk thick- 



ness. The linear analysis is confirmed with detailed simu- 
lations of disk instability, which also show that a similar 
criterion applies to the discussion of non - axisymmetric 
pertu rbations (see, e.g., the review by iDurisen et al.l 
2007). In computing Q we identify the velocity dispersion 
c s with the adiabatic sound speed and we set n = vk/t 
appropriate for Keplerian motion. Both in simulation 
Z4 and in simulation Z4NOMOL the Toomre parameter 
Q > 1 for all radii e < r < 0.3 kpc that cover the two 
disks. In both simulations the inner disk is characterized 
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by Q < 2 while the outer disk is characterized by Q > 2. 

The measured Toomre Q parameters imply that the 
disks in simulation Z4 are somewhat less stable than the 
disks in simulation Z4NOMOL. This is consistent with 
the observation that spiral arms in Z4 are more distinct 
than in Z4NOMOL (right panels in Figures [3] and |4]) and 
is likely a direct consequence of the fact that in simu- 
lation Z4 the disk temperatures are significantly lower 
than in simulation Z4NOMOL thanks to efficient low- 
temperature cooling by molecular hydrogen. Note that 
the stability of the outer and of the inner disk in simula- 
tion Z4 and the stabilitiy of the inner disk in simulation 
Z4NOMOL may be artificially increased due to the im- 
posed Jeans floor as mentioned above. Note also that in 
equating the velocity dispersion with the sound speed, 
we may have underestimated the true velocity dispersion 
and hence the stability of the disks. 

4. DETECTING THE FIRST GALAXIES WITH JWST 

One of the main science goals of the upcoming 
JWST is the detecti on of light from the first galaxies 
(|Gardner et al.ll2QQ6h . Here we present an estimate of 
the expected flux from the first galaxies based on our 
simulations of galaxies with halo masses M v i r ~ 10 9 M 
at z > 10 and investigate their detect ability with the 
instruments aboard JWST. 

JWST will observe the first galaxies using deep field 
imaging and spectroscopy with the Near Infrared Cam- 
era (NIRCam) and the Mid Infrared Instrument (MIRI) 
and using spectroscopy with the Near Infrared Spec- 
trograph (NIRSpec) and als o MIRI (for an overv iew 
of th ese instruments see http://www.stsci.edu/jwst; see 
also iGardner et aLll2006l ). NIRCam will allow imaging 
and low resolution (R = A/ A A < 100) spectroscopy 
within a field of view of 2.2' x 4 A' and an angular reso- 
lution of ~ 0.03" — 0.06" in the range of observed wave- 
lengths A = 0.6 — 5 fim. The multi-object spectrograph 
NIRSpec will enable medium resolution (R ~ 100 — 3000) 
spectroscopy of up to ~ 100 objects simultaneously 
within a field of view of 3.4' x 3.4'. NIRSpec will operate 
in the same wavelength range as NIRCam but at lower 
angular resolution (~ 0.1"). Finally, MIRI will comple- 
ment NIRCam and NIRSpec by providing imaging, low 
and medium resolution spectroscopy within the range of 
observed wavelengths A Q = 5 — 28.8 fim and fields of 
view and angular resolutions of, respectively, ~ 2' x 2' 
and -0.1" -0.6". 

4.1. Intrinsic Luminosities 

Our estimates of the observability of the first galax- 
ies are based on our simulations of galaxies in halos 
with masses ~ 10 9 M at z > 10. We examine and 
compare two idealized scenarios for star formation de- 
rived from the gas accretion rates observed in our sim- 
ulations, chosen to bracket the range of likely scenar- 
ios and parametrized such as to enable the straightfor- 
ward rescaling and extrapolation of our results. We com- 
bine assumptions about the nature of the forming stel- 
lar populations with population synthesis models to es- 
timate the luminosities in the Lya, Ha and Hell (rest- 
frame wavelength A e = 1640 A; hereafter Hel640) neb- 
ular recombination lines and the intensities of the non- 
ionizing (combined stellar and nebular) UV continuum 



(rest-frame wavelength A e = 1500 A; hereafter UV1500). 
We translate line luminosities and UV continuum inten- 
sities into observed fluxes and compare them with the ex- 
pected flux limits for observations with JWST. Based on 
extrapolation of our results to galaxies with both lower 
and larger halo masses, we estimate the number of high- 
redshift halos JWST will detect. 

The first of the two star formation scenarios explored 
here assumes that stars form in a single central instan- 
taneous burst with total stellar mass 

m. = io« m (A) (_^_) , (2) 

where /cool is a conversion factor that determines the 
amount of gas mass available for starbursts inside ha- 
los with virial masses M v i r , and f* is the star forma- 
tion efficiency, i.e., the fraction of the available gas mass 
that is turned into stars. Setting f coo \ = 0.01, this sce- 
nario is motivated by the rapid accretion of large gas 
masses (M g > 10 7 Mq) onto the central unresolved 
core observed in our simulations of halos with virial 
masses M v i r ~ 10 9 M at around z < 12 (see the 
right panel of Figure [2j). We choose a relatively high 
star formation efficiency, /★ = 0.1 , expected for the in i- 
tial bursts (e.g, iWise fc Cer] 120091 : Uohnson et aLll2009D . 
The adopted conversion factor f coo \ = 0.01 between the 
gas mass available for star formation and the halo virial 
mass is consistent with gas collapse fra ctions in previous 
simulations of the first g alaxies (e.g., IWise et all 120081 : 
iRegan fc Haehneltll2009af ). 

The second scenario for star formation considered here 
assumes that stars form continuously at a rate 

M^) = O.O5M 0yr -'(A.)(_|^ (3) 

in proportion to the rate M g at which gas is accreted. 
Indeed, galaxies with masses > 10 9 M may be suf- 
ficiently massive to sustain a moderate level of near- 
continuous star for mation despite ongoing feedback (e.g., 
Wise & Cen 2009; we discuss the effects of feedback in 
more detail in Section [6] below). We adopt a gas ac- 
cretion rate M g ~ 1 M yr _1 that describes the rate 
of accretion of gas onto the central region with radius 
r < 0.1r v i r for z < 15 in our simulations of halos with 
virial masses ~ 10 9 M (see the right panel of Fig. [6]). 
At this radius the gas surface density is roughly in agree- 
ment with the critical surface density > 10—100 M pc -2 
for star formation in the low-redshift universe (see the 
middle panel of Figure [9} ; this critical surface density is 
further discussed in Section [5] below. We approximately 
include the effects of feedback from star formation by em- 
ploying a lower star formation efficiency f* = 0.05 than 
used for the starburst. The implied star formation rates 
M*(z) = 0.05 M yr _1 are consistent with star forma- 
tion rates found in recent feedback simulati ons of galax- 
ies inside halos with masses ~ 10 9 M (e. g. , IWise fc Cenl 
2009; iRazoumov fc Sommer-Larsenll2010h . 

To compute the luminosities of the stellar populations 
that form in the two scenarios we must specify the metal- 
licities of the stars, the stellar ages, and the distribution 
of stellar masses at the time of formation, i.e., the initial 
mass function (IMF). We first assume that starbursts oc- 
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Fig. 10. — Luminosities in the Ha (left panel, left axis) and Lya (left panel, right axis) and Hel640 (middle panel) nebular recombination 
lines and the intensity of the combined stellar and nebular UV1500 continuum (right panel) for galaxies inside halos with mass M V - 1Y = 
10 9 Mq. The estimates are based on the Schaerer ( 2003) stellar population synthesis models and assume Case B recombination theory and 
zero escape fractions. Red diamonds assume the formation of metal- free stars with top-heavy IMF in an instantaneous burst of total stellar 
mass M* = 10 6 Mq at z = 12. We also show the luminosities expected for identical bursts but assuming metal-free stars with a normal 
IMF (blue triangles) and stars with metallicities Z = 5 x 10 -4 Zq and normal IMF (black crosses). The lines assume the continuous 
formation of stars with metallicities Z = 5 x 10 -4 Zq and normal IMF at a rate 0.05 Mq yr -1 . All results scale linearly with the star 
formation efficiencies of /* =0.1 and 0.05 assumed for, respectively, the starburst (equation [2]) and the continuous star formation scenario 
(equation [3])- Note the large dependence of the Hel640 flux on metallicity and IMF. The Hel640 line luminosities estimated for the 
continuous star formation scenario are not shown because they are too low to fall inside the plot range. 



cur in metal-free gas and form clusters of zero-metallicity 
stars. We adopt a top-heavy IMF, i.e., an IMF biased 
towards massive (M* ~ 100 Mq) stars. Such an IMF 
is expected to characterize the first, metal- free genera- 
tion of stars which form by radiative cooling from colli- 
sionally excited m olecular hydrogen (e.g., the review by 
iBromm et aTll2009D . 

The IMF of metal- free stars, however, is still subject 
to large theoretical uncertainty. Stars forming out of gas 
with an elevated electron fraction, such as produced be- 
hind structure formation or supernova (SN) shocks or as 
inside and near ionized regions, could have characteristic 
masses substantially less than ~ 100 M . This is because 
the increased electron abundance boosts the production 
of HD which enables gas to cool to much lower temper- 
ature s than is possible with molecular hydrogen alone 
(e.g., Nakamura & Umemura 2002; Nagakura & Omukai 
2005; Johns on fc Brom m 2006; Stacv & Bromm 20071 : 
see also lShapiro fc Kandll987f and lClark et al.l2010t ). We 
therefore repeat our analysis assuming the formation of 
metal-free stars with a normal IMF, similar to the one 
used to describe star formation in the nearby universe. 

Enrichment to critical metallicities as low as Z c < 
10~ 6 - IO -3 5 Z Q , where we set Zq = 0.02, will 
also imply the tra nsition from a top-heavy IMF to a 
norm a l IMF (e.g.. [Bromm et all 120011: IBromm fc Loebl 
l200l ISchneider etHI l200l ISmith et al.l 120091) . We 
therefore complement our study of metal-free starbursts 
with the study of a burst of stars with above-critical 
but low metallicities Z > 3 x 10" 4 Zq and nor- 
mal IMF. Note that a few massive star SN explosions 
may already be sufficient to enrich the first galaxies 
to metallicities Z > Z c (e.g.. iScannapieco et "all [ 2003: 
jTornatore eF al. 2007; Wise & Abel 2008; Karlsson et all 
120081 : iGreif et al.l 120101 : iMaio et al.l I2010D . We there- 
fore always adopt above-critical (but low) metallicities 
Z > 3 x 10 -4 Zq and a normal IMF in the continuous 
star formation scenario. 

We compute the stellar ionizing luminosities expected 



for the two star formation scenarios using the population 
synthesis models for zero-age i nstantaneous bu rsts and 
continuous star formation from ISchaererl (j2003l, some of 
these models have been previ ously published in ISchaererl 
2002). The Sch aererl (j2003f ) models assume a power- 
law IMF p(ra*) oc m~ a with iSalpeterl (|1955h exponent 
a = 2.35 but allow for different r anges for the m asses 
m+ of individual stars. We use the ISchaererl (2003) zero 
metallicity models for instantaneous starbursts with ini- 
tial masses in the range 50 — 500 Mq and 1 — 100 M 
to describe metal-free stars with a top-heavy and normal 
IMF, respectively. We des cribe the popul ations of low- 
metallicity stars using the ISchaererl j2003) models with 
initial masses in the range 1 — 10 M^ and meta llici- 
ties 6 Z = 5 x 10~ 4 Zq. Following ISchaererl (12003ft. we 
use case B recombination theory (e.g., Ost erbrockl ll989) 
to relate the ionizing luminosities of the stellar popula- 
tions of specific age, mass, and metallicity to the nebular 
luminosities of the surrounding gas, assuming that all 
ionizing photons are absorbed, i.e., that the fraction of 
ionizing photons escaping the galaxy, / esc , is zero. 7 We 
will discuss these assumptions at the end of this section. 

Figure [10] shows the luminosities of the Ka and Lya 
(left panel) and Hel640 (middle panel) nebular lines and 
the intensities of the non-ionizing stellar and nebular 
UV1500 continuum (right panel) expected for the two 
star formation scenarios described above. Note that the 
luminosities in the Lya line are related to those in the Ha 
line by t he simple scali ng L(Lya) w 8.6L(Ha) (Tables 1 
and 4 in Schaerer 2003). For low metallicity and normal 
IMF, the starburst scenario implies line luminosities and 
continuum intensities that are roughly twice as large as 

6 Our conclusions are insensitive to the precise choice for the 
metallicity Z > in the models. 

7 We obtain the luminosities in the nebular recombination lines 
using equations (7) and (8) in Schaerer (2003) together with the 
data in their Tables 1, 3, and 4. We obtain the combined stellar 
and nebular UV1500 continuum intensity, averaged within a 20 A 
band center ed on 1500 A, from the corre sponding online data sets 
provided at http://obswww.unige.ch/sfr 
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those implied by the continuous star formation scenario. 

At fixed IMF, the luminosities in Ha and Lya and the 
UV1500 continuum intensities of the starbursts are in- 
sensitive to the stellar metallicity. The zero-metallicity 
starburst models with top-heavy IMF imply Ha and 
Lya line luminosities and U VI 500 continuum intensities 
larger by about one order of magnitude than those im- 
plied by the zero-metallicity starburst model with a nor- 
mal IMF. In contrast, the starburst luminosities in the 
Hel640 line depend strongly on both the IMF and the 
stellar metallicity. At fixed normal IMF, a change from 
low to zero stellar metallicity implies an increase in the 
He 1640 line luminosity by about three orders of magni- 
tude. This is because the exceptionally hot atmospheres 
of zero-metallicity stars turn them into strong emit- 
ters of Hell ionizing rad i ation fe.g., ITumlinson fc Shulll 
I2QQQI : iBromm et al.ll200TI : ISchaererll2003f ). An additional 
change from normal to top-heavy IMF increases the lumi- 
nosity in the He 1640 line by another order of magnitude. 

The large differences in He 1640 line luminosities of- 
fer the prospect of distinguishing observationally be- 
tween stellar populations made of metal-free and metal- 
enriched stars and of constraining their IMFs (e.g., 
Tumlinson fc Shulll [2QQQt IBromm et al.l 12001 IQh et all 



20011 : Uohnson et all l2QQ9f T The Hel640 recombina- 
tion line will also be excited due to the emission of 
ionizing radiatio n from a central accreting black hole , 
if present (e.g., IQh^taD 120011 : ITumlinson et al.l 120011 : 
Uohnson et al.l |201Q[ ) . Observationally, accreting black 
holes could be distinguished from metal-free stellar popu- 
lations through the de tection of their X-ray emission (e.g, 
lHaiman fc Lo eb 1999). Note though that an evolved stel- 
lar populat ion may also contribute to the X-ray emis- 
sivity (e.g.. IOl3l200TI : iPower et ^ 12009). X -ray sources 
may ionize the gas in a larger region than stellar sources, 
implying a spatially more extended emission of recombi- 
nation radiation. 

4.2. Observed Fluxes 

We translate the line luminosities and UV continuum 
intensities into observed fluxes to investigate the de- 
tectability with JWST. The flux density from a spatially 
unresolved object emitted in a spectrally unresolved line 
with rest- fram e wavelength A e and line lum inosity L is 
by (e.g., [Qhl[T999L ilohnson et aHl2009h 



given 



/,(A ) = 



LX e (l + z)R 
4ircdl(z) 



(4) 



' 3 nJy 



L 



10 40 erg s- 1 



/ A e \ / R \ n+z 

V1216 Ay Vioooy V 11 



where A e is the rest-frame wavelength, A Q = (1 + z)X e 
the observed wavelength, and d\, ~ 100 [(1 + 2)/ 10] Gpc 
the luminosity distance. The UV continuum intensity L v 
of a spatially unr esolved ob ject implies an observed flux 
density (e.g.. IOhllT99l IBromm et al.ll200lD 
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Fig. 11. — Stellar masses M* jm i n of the lowest mass starburst 
observable through the detection of the Ha line (solid curves) or 
the Hel640 line (dashed curves) or the UV1500 continuum (dash- 
dotted curves) with JWST, assuming an exposure of t ex p = 10 6 s 



10. The masses scale with f ( 



-1/2 



Stellar masses derived 



andS/N- xu. xu 
from the ISchaererl (2003) zero-metallicity starbursts with normal 
IMF, the zero-metallicity starbursts with top-heavy IMF, and the 
low-metallicity starburst, are shown, respectively, in blue, red, and 
black. The right axis shows the masses M m [ n = 10 3 M* >m i n of 
halos expected to host a starburst with stellar mass M* )m j n . The 
conversion between stellar and halo masses is based on equation ([2} 
with /cool — 0-01 an d /* = 0.1. For reference, the dotted curve 
shows the virial mass (with corresponding labels on the right axis) 
for a halo with virial temperature T v i r = 10 4 K. The black dashed 
curve is not shown because it exceeds the plot range. 

Flu x densities , f v , are related to AB magnitudes, ttiab, 
via (|Oki[l974; O ke fc GunnlflQSl 



rriAB 



-2.5 log 



10 



nJy 



31.4. 



(6) 



The assumption that the lines are spectrally unre- 
solved is excellent for both Ha and He 1640, whose line 
widths AA/A < 10- 4 (T/10 4 K) 1 / 2 are set by thermal 
Dopp ler broadening at temperature T < 10 4 K (e.g.. lOhl 
1999). We also note that at redshifts z > 10 a trans- 
verse physical scale Al corresponds to an observed an- 
gle A0 = Al/d A - 0.1"(AZ/0.5 kpc)[(l + z)/10], where 
dx = (1 + z)~ 2 d\, is the angular diameter distance. 
Hence, if most of the nebular emission originates from 
within the vicinity of the stellar populations, which we 
assumed to be concentrated in the halo centers, i.e., at 
r < 0.1r v i r , the assumption that the emitting regions 
are spatially unresolved is good for both the Ha and the 
He 1640 line and it applies equally well to the UV contin- 
uum. 

In contrast, the Lya line radiation undergoes resonant 
scattering. Henc e, it will likely be additionally spectrally 
broadened (e.g., iNeuf eld 1990), and sp atially extended 
with typical angular size AO - 15" ([Loeb fc Rvbickil 
Il999[ ) . It will be damped d ue to absorption by intergalac- 
tic ne utral hydrogen (e.g.. iMiralda-Escudel 1998; S antosl 
2004). Note that Lya radiation from galaxies at red- 
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Fig. 12. — Observed flux densities in the Ha and Hel640 recombination lines and of the combined stellar and nebular UV1500 continuum 
derived using the line luminosities and continuum intensities shown in Figure [Tol and using equations @J and (|5]). The flux densities scale 
linearly with the star formation efficiencies of /* = 0.05 and 0.1 for, respectively, the continuous star formation and the starburst scenario. 
Dotted lines show the sensitivity limits for observations with JWST, assuming exposures of 10 4 , 10 5 , and 10 6 s (top to bottom) and 
S/N=10. The Hel640 line fluxes obtained for the continuous star formation scenario are not shown because they are too low to fall inside 
the plot range. With exposures texp ^ 

10 6 s, JWST will have the sensitivity to distinguish between metal-free starbursts with top-heavy 
IMF and metal-free or metal-enriched starbursts with normal IMF inside ~ 10 9 M© halos based on the detection of the Hel640 line. 



shifts z > 10 may be particularly strongly damped be- 
cause the reionization of the universe was pro bably only 
accom plished at much lower redshifts (e.g., iFan et all 
2006). On the other hand, scattering off outflowing 
inters tellar gas may help the Lya radiation to escape 
(e.g., iDijkstra fc Wvithel l2010h . and galaxies may re- 
side in a ionized bubble sufficiently large for Lya pho- 
tons to redshift away in the expand i ng universe (e.g. , 
| Cen fc Haimanl [2000t iHaimanl l200l lLoeb et al.l [20051: 
IWvithe fc Loebl 120051 : iLehnert et al.ll2QlQh . Since we do 
not treat radiative transfer effects here, Lya line fluxes 
implied by equation ((4]) must be considered upper limits. 
In the following we therefore mostly discuss the observ- 
ability of the Ha and Hel640 lines and of the UV1500 
continuum. 

The Hel640 recombination line (A e = 1640 A), as well 
as the Lya line (A e = 1216 A) not further discussed here, 
will be detected by JWST with NIRSpec at a spectral res- 
olution R ~ 1000, while the Ha line (A e = 6563 A) will 
be detected with MIRI at a spectral resolution R ~ 3000. 
JWST will detect the UV1500 (A e = 1500 A) continuum 
using NIRCam. As an illustration, Figure [TT] sho ws the 
minimum stellar mass required for the Schaerer (2003) 
model starbursts employed here to be observable with 
JWST. We assume exposures with signal to noise ratio 
S/N = 10 and duration t exp = 10 6 s and the currently ex- 
pe cted flux limits 8 for o bse rvations with J WST (Table 10 
in lGardner et al.ll2006l; seelPa nagia 2005 for a graphical 
presentation and [Johnson et al.l 120091 for a useful sum- 
mary) . Figure [11] demonstrates that even for exposure 
times as long as 10 6 s, JWST wi\\ not have sufficient sen- 
sitivity to detect stellar populations with masses below 
~ 10 5 -10 6 M©. JWST will thus not be a ble to see stellar 
light from individual first stars (e.g.. lOh| [T999: Oh et al.l 
120011 : iGardner et al.ll2006l ). 

Figure fT2l shows the Ha (left panel) and Hel640 (mid- 
dle panel) recombination line fluxes and the non-ionizing 

8 Flux limits /ij m reported in Gardner et al. (2006) assume S/N 



10 and t e 



10 4 s. Here we rescale these limits to other 



exposure times using fn m oc t, 



1/2 

exp 



Flux limits for t e 



10 6 



and S/N=10 can be read from Figures \12\ and ITTTI 



UV1500 continuum fluxes (right panel) for the models 
presented in Figure fT()l The JWST flux limits for expo- 
sure times t exp = 10 4 , 10 5 and 10 6 s are indicated by the 
dotted lines in each panel. The figure reveals that the 
scenario of continuous star formation implies line and 
continuum fluxes too low to be observable, even when 
assuming exposure times as large as 10 6 s. JWST, how- 
ever, may see starbursts similar to those modelled here. 
In exposures with duration ~ 10 6 s, MIRI will detect 
such starbursts in Ha for all metallicities and IMFs ex- 
plored. 

Figure [12] also reveals that JWST has the potential 
to constrain the properties of starbursts in galaxies with 
halo masses as low as ~ 10 9 M©, based on the detec- 
tion of the Hel640 line. Indeed, only the zero-metallicity 
starburst with a top-heavy IMF and observed with an 
exposure of < 10 6 s is detected in Hel640. Starbursts 
inside > 10 — 100 times more massive halos would be de- 
tected in the Hel640 line independent of whether their 
IMFs are top-heavy. Their nature could then be further 
constrained by measuring the ratio of the Ha and Hel640 
line strengths. 

Note, finally, that if scattering by the intergalactic gas 
can be ignored, the total (i.e., integrated over the line) 
flux in the unresolved Lya line would be a factor ^8.6 
larger than the total flux in the unresolved Ha line, based 
on the relation L(Lya) ~ 8.6L(Ha) noted above and 
shown in Figure [TOl For a galaxy at redshift z « 10, 
JWSTs NIRSpec is a factor of « 3 more sensitive to 
the detection of the redshifted Lya line than MIRI is to 
the detection of the redshifted Ha line 9 . That is, unless 
radiative transfer effects cause JWST to see less than 
1/(8.6 x 3) « 4% of the Lya line flux, the Lya line will 
be easier to detect than the Ha line. Hence, despite the 
large uncertainties arising from its resonant nature, the 
Lya line re mains a powerful probe of h igh-redshift galaxy 
formation (jPartridge fc Pe ebles 1963). 

4.3. JWST Number Counts 

9 See the sensitivity limits for detection 

of narrow unresolve_d line fluxes quoted at 

http: / / www.stsci.edu/jwst / science/data_simulation_resources / sensitivity 
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Fig. 13. — JWST starburst counts. Left panel: Number of halos N(> z) with redshifts > z and masses > M m i n , where M m i n is the 
lowest mass halo capable of hosting a starburst observable through the detection of the Ha line (solid curves) or the Hel640 line (dashed 
curves) or the UV1500 continuum (dash-dotted curves) with JWST (Figure [TTJ right axis). We have assumed an exposure of t e xp = 10 6 s 
and S/N = 10. Counts for the model starbursts of zero metallicity and normal IMF, zero metallicity and top-heavy IMF, and low non-zero 
metallicity are shown, respectively, in blue, red, and black. The black dashed curve is not shown because it falls below the plot range. 

1/2 

Right panel: Number of halos N(> /) above z > 10 with observed fluxes > /. The vertical lines show the JWST flux limits fy im oc t e ^p 
for observations of the YLa line (solid), the Hel640 line (dashed) and the UV1500 continuum (dot-dashed), assuming exposures t exp = 10 6 s 
and S/N = 10. The black dashed curve is not shown because it falls below the plot range. JWST may detect a few tens (for nonzero 
metallicities and normal IMFs) up to a thousand (for zero metallicity and top-heavy IMFs) starbursts with redshifts z > 10 in its field of 
view of ~ 10 arcmin 2 . Our estimates for the number counts scale linearly with the assumed starburst durations of r s b = 3 Myr and 30 Myr 
for, respectively, the starbursts with top-heavy and normal IMFs. 



How many star- forming galaxies can we expect JWST 
to detect? For simplicity and brevity of the presentation 
we ignore that galaxies may form stars in a continuous 
mode and assume that starbursts shine at constant lumi- 
nosity over a time interval r s b- The number of galaxies, 
per unit solid angle, above redshift z that JWST will 
detect is then obtained using 

go,,. r*>™ -s*r iu*a^i 

(7) 

where tu(z) is the age of the universe at redshift z, 
dV — cH~ 1 (z)di j (l + z)~ 2 is the comoving volume el- 
ement, H(z) = H [n m (l + zf + ^ A ] 1/2 , and H = 
100k km s _1 Mpc -1 is the Hubble constant. We ap- 
proximate the comoving number density n(M, z) of halos 
with mass M at redshift z by the P r ess- Schechter halo 
abundanc e dPress fe Schechterl [19741 : iBond et al.l 119911 : 
see, e.g., Zentnerl [20071 for a recent review). We set 

M min (z) = f~*oif* lM *,mia( z )i wnere M* )min (» is the 
smallest stellar mass observable at redshift z (Figure [TT]). 
and we use / CO oi = 0.01 and /* = 0.1 (see equation [2]). 

Figure [13] shows our estimates of the number of ob- 
servable starbursts in exposures of t exp = 10 6 s and 
S/N=10. The estimates scale linearly with the assumed 
durations r s b of the starbursts. We set r s b = 3 Myr 
for the zero-metallicity starbursts with top-heavy IMF, 
which approximately corresponds to the time it takes for 
its massive ~ 100 M stars to age and explode, upon 
which further star formation, if not suppressed by SN 



feedback, will more likely occur inside metal-enriched 
gas, hence ceasing the zero metallicity burst. We set 
r s b = 30 Myr for the zero-metallicity and low-metallicity 
starbursts with normal IMF. Our choice for this longer 
duration reflects the longer time it takes, on average, for 
massive stars forming inside bursts with normal IMFs to 
evolve and explode in SNe. Starburst durations up to 
~ 30 Myr are consistent with the gas dynamics and the 
amount of gas available for star formation in our simula- 
tions (see Section 13^21 and Figure [2]). 

Figure [13] demonstrates that JWST will enable the 
detection of a few tens up to a thousand star-bursting 
galaxies with redshifts z > 10 in its field of view of 
~ 10 arcmin 2 . JWST will allow to constrain the pre- 
dominant nature of the first starbursts as the Hel640 
recombination line is only detected in significant num- 
bers for the case of zero-metallicity starbursts with top- 
heavy IMF. Intriguingly, our estimates imply that the 
first galaxies may be more readily detectable in Ha spec- 
troscopic searches than in UV continuum surveys (see 
also Figure [T2]) . The expected number of star-bursting 
galaxies with redshifts z > 10 to be detected with JWST 
in exposures other than 10 6 s can be read from the right 
panel in Figure [13] Our estimates are consistent with 
previous estimates of JWST starburst counts for similar 
assumptions ab out the conversion between st ellar and 
halo mass (e.g., lHaiman fc Loeblll99l IOhlll999[ ). 

The estimates of the observability of the first galax- 
ies are uncertain due to the assumptions underlying the 
computa tion of UV lin e and continuum emission. We 
followed Schaerer (2003) and computed the nebular con- 



THE FIRST DISK GALAXIES 



15 



tribution to the stellar luminosities ol the first galaxies 
using Case B recombination theory. iRaiter et al.l (j2010l ) 
point out that the Case B approximation ignores pho- 
toionizations from excited states and hence underesti- 
mates the photoionization rate. At low nebular metallic- 
ities and for hot stellar sources, their detailed photoion- 
ization models imply Lya line luminosities and nebu- 
lar UV continuum intensities larger by factors of 2 — 3 
(see their Fi gure 10) than expe cted under the Case B 
assumption. IRaiter et al.l (|2010f ) also find that the line 
luminosities in Ha are insensitive to whether Case B is 
assumed. Hence, Lya line and UV continuum emission 
may provide relatively stronger signatures of star forma- 
tion inside the first galaxies than suggested here. 

Finally, we have assumed that a negligible fraction / esc 
of stellar ionizing photons escapes the star-forming re- 
gions without being converted into recombination radia- 
tion by the surrounding gas, i.e., f esc = 0- Recent numer- 
ical work has emphasized that the escape fraction may 
depend strongly on the structural properties of galax- 
ies as determined by their masses and internal processes 
like star formation and feedback (e . g., Fuiita et al. 2003; 
Gnedin et all l2008ri IWise fc Cen 2009;" Johnson et all 
2009; iRazoumov fc Sommer-Larseri 1201(1 lYajima et al.l 
2010). While escape fractions of order unity f esc ~ 0.5 
are possible for low-mass (< 10 9 M ) halos with turbu- 
lent gas dynamics and amorphous morphology, signifi- 
cantly smaller escape fractions are expected for galaxies 
that form most of th eir stars inside dense rotationally- 
supp o rted disks (e.g., iGnedin et al1l2008ri IWise fc Cenl 
2009; IRazoumov fc Sommer-Larsenl I2010D . The lumi- 
nosities implied by our assumption of a zero escape frac- 
tion can be rescaled to account for non-zero escape frac- 
tions by multiplication with the factor (1 — / esc )- 

5. DISCUSSION 

An interesting outcome of our simulations is the col- 
lapse of the halo gas into two extended rotationally 
supported disks. When and how the first galaxy- 
scale disks were formed is currently not well under- 
stood. Extended disks are found in large-scale hy- 
drodynamical cosmological simulations and in cosmo- 
logical simulations of individual massi ve (> 10 1Q M ) 
halos o r halos at low red s hifts (e.g., Kaufm ann et al.l 
20071: iMashchenko et~aT1 120081: Levine et all 120081: 
Sawala et al. 2010; G overnato et al.l 120101 : iSchave et al.l 



20101 : iSales et al.ll2010l ). Cosmological simulations of the 
first minihalos and low-mass (< 10 8 M^) halos, how- 
ever, have not yi e lded such disks (e. g., Wise et al. 2008; 
Greif et a l. 2008; Wi se fc Cenll200l [Regan & Haehiitl 



2009al: iStacv e t al. 2010|). The gas dynamics inside these 
halos is dominated by turbulent motions instead. This 
morphological bimodality suggests that mass is an im- 
portant factor in determining whether a given halo may 
host an extended disk. Evidence for the suppression of 
disk formation in low-mass galaxies comes from observa- 
tions of dwarf galaxies in the local universe which suggest 
a critical stellar mass belo w which stellar disks become 
systematically thicker (e.g ., iSanchez-Janssen eTaH 1201(1 
iRovchowdhurv et al.ll20Tof ). 

Our simulations show that the formation of extended 
gas disks is possible in halos with masses as low as 
M vi r ~ 1Q 9 Mp at redshifts as high as z > 10 (but 
see lLatif et al.l l20To[ ). However, we acknowledge that 



the assembly of the disks may be determined in part 
by the imposed Jeans floor. The Jeans floor artificially 
heats the gas and increases the sound speed, and this re- 
duces the Mach numbers at which accretion shocks sup- 
ply fresh gas to the central high-density regions. Our 
simulations may therefore potentially underestimate the 
ability of accretion flows to stir up the gas and to chan- 
nel energy into turbulent motions, which could other- 
wise impede the formation of thin extended disks by 
transporting angular momentum and driving material 
that would without turbulence have circularized in a thin 
disk rapidly inward. Indeed, high Mach number accre- 
tion flows have been identified as major drivers of the 
turbulent gas motions and morphologies seen in simu- 
lations of high-redshift low-mass galaxie s (jWise fc Abell 
120071 : IWise et al.ll2008l : IGreif et al.ll2008h . Note that in 
the simulation without molecular cooling (ZJ±NOMOV), 
the pressure floor affects only the inner disk. The for- 
mation of the outer disk thus is a robust outcome of this 
simulation. 

The unresolved gas cores embedded within the disks in 
our simulations have masses (~ 5 x 10 7 M , see Figure |2j) 
and compact sizes (< 10 pc) to potentially evolve into nu- 
clear star clusters or massive black holes in the centers 
of dwarf galaxies. There is an observed relation between 
the masses M ga i of galaxies and the masses M cent of the 
nuclear star clusters in spheroidal galaxies or of mas- 
sive black holes in ellipticals and bulges, both of which 
are M cent ~ 2 x 10~ 3 M ga1 (e.g., iFerrarese et al.l 12001 
IWehner fc Ha rris 2006|, and references therein). It is not 
straightforward to define the masses of the galaxies in our 
simulations in a manner consistent with the observational 
definitions (see, e.g., the discussion in Li et al. 2007). 
However, we may conservatively assume that galaxy 
masses M ga i are fractions < 1 of the associated halo 
virial masses M v i r . Then, the masses of the central cores 
are too large by factors > (M cen t/M vir )/(2 x 10 -3 ) = 
(5 x 10 7 M ? /10 9 M )/(2 x 10- 3 ) = 25 to fit the ob- 
served relation. Feedback from a central source, which 
was ignored here, may be crucial for establishing this 
relati on (e.g., McLaughlin et al. 2006; Naravanan et al.l 
2008; Jo hnson et al.l 120101 : but see, e.g.. iLi e 
Larson! I2010D . 
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However, it is not known if the locally 
observed relation is already established at the high red- 
shifts of interest and if it extends to the low-mass halo 
regime considered here. 

We can combine the surface density profiles of the disks 
obtained in our simulations with assumptions about a 
threshold surface density for star formation to speculate 
on the stellar radii of the simulated galaxies. A thresh- 
old of ~ 10 M pc -2 is implied by observations at kilo- 
parse c scales of star formation in nearby disk galaxies 
(e.g lKennicuttlll989l : iKennicutt et al]ll998t iBigiel et al.l 
2008) and is s upported by semi-analytical and numeri- 
cal work (e.g.. lElm e green & Parravano 1994; E lmegreenl 
120021 : ISchavdl2003 IGnedin fc Kravtsovl l20Toh . At high 
redshifts this threshold could be larger, < 100 M pc -2 , 
mostly because of the low dust abundances (implying less 
shielding from the supposed UV background) a t these 
epochs (|Gnedin fc Kravts ov 2010; see also, e.g., iSchavd 
2004; Krumh olz et al.ll2QQ9h . Our simulations then imply 
stellar radii < 0.1 kpc (see Figure |9j). Such small stellar 
radii are characteristic of dwarf-globular transition ob- 
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jects and small dwarf spheroidals around the Milky Way 
a nd other members of the Local Group (see Figure 8 
in lBelokurov et a"T1l2007[ ). We, however, caution that the 
relatively massive dwarf galaxies simulated here may con- 
tinue their stellar growth well below z < 10 as they may 
accrete gas also after reionization has raised the Jeans 
mass in the intergalactic medium to ~ 10 8 M (see Sec- 
tion [6] below) . The possibility remains that the central 
gas core forms stars but the disks do not, in which case 
the galaxies would, upon gas loss, evolve into a massive, 
compact star cluster. 

We have shown that the detection of recombination 
lines and UV continuum radiation emitted by gas sur- 
rounding the first stellar populations in deep exposures 
with the upcoming JWST will likely only probe galax- 
ies inside halos with masses > 10 9 M©. Detection of 
stellar light and recombination radiation from smaller 
galaxies may be pos sible if these galax ies are gravita- 
tionally lensed (e.g., Uohnson et al]|2009[ ). Note that re- 
combination radiation may also be produced by halo gas 
that does not join the central disks smoothly but comes 
to a halt in a shock. The infall energy of the shocked 
gas would be transformed into radiation, ionize the disk 
environment and be re-e mitted as recombination lines 
(jBirnboim fc Dekel 2003). We have ignored this poten- 
tially significant contribution to the recombination line 
luminosities in our study of the observability of the first 
galaxies presented here. 

In addition to detecting recombination radiation from 
the interstellar gas, JWST may observe the first galaxies 
through the detection of cooling radiation emitted during 
their assembly. JWST will probably not have the sensi- 
tivity to detect Lya cooling radiation fro m galaxies with 
halo virial masses M v j r < 10 10 M (e.g., Hai man et al.l 
120001 : lDiikstra![2009h . However, JWST and other tele- 
scopes such as the Atacama Large Millimeter Array or 
the proposed Single Aparture Far-Infrared Observatory 10 
and Space Infrared Telescope for Cosmology and Astro- 
physics 11 may detect these low-mass galaxies through 
the cooling radiation e mitted by molecular hydrogen (for 
an overview see, e.g., lAppleton et al.ll2009f ). The lumi- 
nosity in cooling radiation from molecular hydrogen will 
be strongly boosted if emitted by gas inside SN shells 
(jCiardi fc Ferrarall2001f ) or by gas powered by X-ray ir- 
radia tion from a central black hole (Spaans & Mei ierinkl 
l2008h . 

6. LIMITATIONS AND FUTURE WORK 

In this work we presented our first steps towards self- 
consistent simulations of the formation and evolution of 
the first galaxies. As such, we have limited ourselves to 
the study of important aspects of the gravitational as- 
sembly of individual dwarf galaxy halos and of the gas- 
dynamical processes inside their virial regions. Our goal 
is to build, step by step, ever more realistic simulations 
that will allow us to draw an increasingly detailed picture 
of high-redshift dwarf galaxy formation. The most im- 
portant challenges for future work concern effects related 
to the formation of stars and associated feedback that we 
have ignored here. Processes that are known to strongly 
affect the assembly and evolution of galaxies include pho- 

10 http://safir.jpl.nasa.gov 

11 http://www.ir.isas.jaxa.jp/SPICA/ 



todissociation, photoionization, SN explosions and asso- 
ciated chemical enrichment and r adiation pressure from 
stellar clusters or black holes (see lCiardi fc Ferr ara 2005 
for a review). 

Our metal-free atomic cooling simulation is consis- 
tent with scenarios outlined in previous works in which 
molecular hydrogen formation and star formation and 
feedback are suppressed in progenitors of the assem- 
bling dwarf galaxy due to the presence o f a photodisso- 
ciatin g Lyman- Werner back g round (e.g., [O h & H aimanl 
[20021; [Johnson et al.l 120081 : [Regan & HaehneH l2009al : 
iRegan fc Haehndtl [20091 IShang et al.l 12010ft . Star 
formation in the progenitors may also be efficiently 
suppressed due to photoionization from early (local) 
reionization, which affects the gas fractions in low- 
mass halos primarily by boiling t he gas out of the 
shallow halo potential wells (e .g.. [Thoul fc Weinberg! 
1996; Barkana fc Loebl [1991 iKitavama et al.l [2000; 
Gnedin 2000; Diikstra et al. 2004; Shapiro et al. 2004; 
IWhalen et aLll2QQ8at iPetkova fc Spring e! 2010). Reion- 
ization also raises the cosmologica l Jeans mass in 
the ionized intergala c tic gas (e.g., [Shapiro et all 119941: 
iGnedin fc Hull pM [Gnedin 2000; Hoeft et aD 120061 : 
lOkamoto etall 120081 : IPetkova fc Springe j l2Ql(fl ) . and it 
affects the rate at which gas can cool a nd sink towards 
the h alo centers (jEfstathiou 1992|; see also|W iersma et al.l 
2009). Both effects lower the star formation efficiency be- 
cause they prevent or impede the replenishing of the gas 
in photoevaporated low-mass halos as well as the accre- 
tion of gas by subsequent low-mass halo generations. 

The effects of photodissociation and photoionization 
on the final state of our simulated galaxies are diffi- 
cult to assess without detailed radiative transfer simu- 
lations, also becaus e internal radiation sources may play 
an important role ([Miralda-Escudel [20051 : iSchavd 2006: 
Gnedin 2010). Preliminary numerical experiments based 
on simulations that include the effects of photodissocia- 
tion and photoionization from a UV background in the 
optically thin approximation 12 indicate that UV radia- 
tion is unlikely to prevent the formation of disks in our 
simulations. However, photoionization may affect the 
disk structure, e.g., because it may determine the lo- 
cal star formation efficiency through its effects on the 
Jeans mass and by producing free electrons that catalyze 
the formation of molecular hydrogen. This latter effect 
could partially offset the negative f eedback on star for- 
mation from photodissociation (e.g . , lHaiman et al.l 2000; 
iRicotti et all 120021 : 1 Ahn fe Shapiro! 120071 ) . 

The assembly and structure of the disks in our sim- 
ulations would probably have been rather different had 
SN explosions been taken into account. Material ejected 
by SNe could sweep up and shock-heat the surround- 
ing gas and entrain strong outflows, even in those rela- 
tively massive halos whose evolution is hardly a f fected 

by photoionizatio n (e. g., [M ac Low & Ferrara 1999; 

Ferrara & Tolstov 2000; Mori et al. 2002; Tassis et al. 
120031 : iDalla Vecchia fc Sch ave 2008|). SNe may drive tur- 

12 We have performed a simulation identical to Z4 but assuming 
equilibrium cooling i n the presence o f a uniform photodissociating 
and photoionizing Haardt & Madau (2001) background in the op- 
tically thin approximation (and with UV inte nsities equal to their 
z = 9 values for all z > 9; see Pawlik et al. 2009 for a detailed 
description of the simulation technique) and using only 3 instead 
of 4 zoom levels, corresponding to 8 times lower mass resolution. 
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bulence and establish a multi-phase interstellar medium 
with hot shock-heated chimneys and c old molecular spots 
inside a warm photoionize d gas f e. g..lMcKee fc Ostrikerl 
19771: IWada fc N orman J^lt^jRicotti et al.l 12008: 
Wise fc AbeJ [20081 : iGreif et al.l [201(1 . The dynamics, 
morphology and thickness of gas disks may then criti- 
cally depend on the distribution of mass over these three 
phases (e.g., Kaufma nn et al.ll2007[ ). Note that feedback 
from SNe could be significantly amplified by previous 
episodes of photoionization (e.g., Kitavama fc Yoshid"al 
120051 : iPawlik fc Schavd 120091 : lHambrick et al.ll2010D . 

SN explosions affect the subsequent star for- 
mation process also by enriching the interstellar 
and intergalactic gas with th e met als synthesized 
in s t ars (e.g., lAguirre et al.l 120011: iMadau et al.l 
[20011: ISclmnameco et al.l 120021 : ICen fc Chisaril 120101 : 
IWiersma et al.l l2Q10h . The increased metallicity en- 
ables additio nal cooling which may help the gas to frag- 
ment (e.g., iBromm e t al. 2001; Schneid eretaTI 120061 : 
ISafranek-Shrader et al.l I2010D . Uappsen et all ([2009D 
have compared high-redshift low-mass halo simulations 
that include low-temperature cooling by both metals 
and molecular hydrogen with identical simulations that 
include only cooling by molecular hydrogen and found 
roughly equivalent levels of fragmentation in both types 
of simulations. Our molecular cooling simulation hence 
may have already captured some of the most important 
effects of low-temperature metal cooling. 

7. SUMMARY 

Motivated by the exciting prospect of the direct de- 
tection of stellar light from redshifts z > 10 with the 
upcoming JWS T, we investigated the assembly of the 
first dwarf galaxies using high-resolution cosmological 
zoomed smoothed particle hydrodynamics simulations of 
individual halos. Previous works suggest that galaxies 
inside halos with masses M v i r < 10 8 M at z > 10 are 
likely too faint, by at least a factor of 10, to be ob- 
served in the proposed exposures with JWST. Hence, 
the light collected in future observations with JWST 
may come mostly from galaxies inside halos with masses 
M vir - 10 9 M . 

We performed two simulations of such galaxies that 
were identical except for differences in the employed non- 
equilibrium primordial gas chemistry and cooling net- 
work. In the first of these simulations, gas cooled by 
emission of radiation from both atomic hydrogen and 
helium and molecular hydrogen. We compared this sim- 
ulation to one in which the formation of molecular hydro- 
gen was suppressed and, hence, the gas cooled only via 
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